Effect of macromolecular crowding on the rate of diffusion-limited enzymatic reaction 
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The cytoplasm of a living cell is crowded with several macromolecules of different shapes and sizes. 
Molecular diffusion in such a medium becomes anomalous due to the presence of macromolecules 
and diffusivity is expected to decrease with increase in macromolecular crowding. Moreover, many 
cellular processes are dependent on molecular diffusion in the cell cytosol. The enzymatic reaction 
rate has been shown to be affected by the presence of such macromolecules. A simple numerical 
model is proposed here based on percolation and diffusion in disordered systems to study the effect 
of macromolecular crowding on the enzymatic reaction rates. The model explains qualitatively some 
of the experimental observations. 
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INTRODUCTION 

The aqueous phase of cell cytoplasm is crowded with 
macromolecules such as soluble proteins, nucleic acids 
and membranes The influence of such crowding 

on biochemical reactions inside physiological media are 
manifold Q. Due to crowding, the average free en- 
ergy fi of a nonspecific interaction between any molecule 
in the medium and a crowding molecule may change 
considerably which may influence the reaction activity 
7 = exp(/i//csT), where ks is the Boltzmann constant 
and T is the absolute temperature. Steric repulsion is 
the most fundamental of all interactions between macro- 
molecules in solution at finite concentration and as an 
effect of such repulsion the macromolecules occupy a sub- 
stantial volume fraction in the cell interior Signifi- 
cant volume fraction of macromolecules in the medium 
imposes a constraint on introducing any new macro- 
molecule. As a consequence of crowding, macromolecular 
association reactions become increasingly favorable. Be- 
cause of crowding, the molecular diffusion in the medium 
is expected to be anomalous Q. The effect of macro- 
molecular crowding on different kinetic steps of enzyme 
catalysis such as formation of enzyme-substrate complex 
and enzyme-product complex were analyzed through dif- 
ferent equilibrium thermodynamic models Q. A number 
of approaches have been proposed to study the effects 
of macromolecular crowding on the reaction kinetic rate 
laws such as the law of mass action fractal like ki- 
netics @, the power law approximation stochastic 
simulation^ and lattice gas simulation^. In these ana- 
lytic and numerical models, the influence of macromolec- 
ular crowding on both equilibrium thermodynamics and 
reaction rates were addressed and it was observed that 
the rate decays exponentially with time as expected in 
equilibrium kinetics. The influence of macromolecular 
crowding on the enzymatic reaction rates has been inves- 
tigated experimentally using a variety of crowding agents 
|lfj |. These studies have also indicated a significant in- 
fluence of crowding on the rate parameters of the enzy- 



matic reaction. It was found that the effect of crowding 
on the enzymatic reaction may be different depending 
on whether the product formation in the enzyme reac- 
tion is limited by the diffusional encounter of substrate 
and enzyme or the formation of the transition state com- 
plex, an association of enzyme and substrate. Moreover, 
molecular diffusion is known to be the major determi- 
nant of many cellular processes and plays a key role in 
cell metabolism where the encounter of the free substrate 
with an active site of the enzyme is often the rate de- 
termining step. However, how the kinetics of an enzy- 
matic reaction is dependent on the size and concentration 
of the crowding macromolecules is still not fully under- 
stood. The macromolecular crowding till date remains 
under appreciated and neglected aspect of the intracellu- 
lar environment [111 ]. It is hence essential to understand 
the experimental observations from microscopic origin. 

In this paper, an approach based on non-equilibrium 
dynamics of enzymatic reactions in the diffusion limited 
regime is considered. The aim is to understand quali- 
tatively the influence of inert macromolecular crowding 
on the diffusion limited enzymatic reactions governed by 
non-equilibrium thermodynamics. A simple numerical 
model in two dimensions (2d) based on molecular diffu- 
sion in disordered systems coupled with enzymatic reac- 
tion is proposed here. The disordered system is mod- 
eled by percolation clusters (T^. It is predicted that 
the rate of a diffusion-limited enzyme-catalyzed reaction 
will experience a monotonic decrease with increase in the 
fractional volume occupancy of the crowding agent. The 
model explains qualitatively certain experimental obser- 
vations. 



THE MODEL 

In brief, the enzyme kinetic reaction in the cell cyto- 
plasm can be described as substrate molecules diffusing 
through crowding macromolecules and binding to the ac- 
tive site of the freely floating enzymes. Subsequently a 
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product is formed if the reaction is energetically favor- 
able and this product diffuses through the same crowd of 
macromolecules. The classical Michaelis-Menten equilib- 
rium enzyme kinetic reaction is given as (l3j 



E + S^ES 



E 



(1) 



where E represents enzyme, S represents substrate, P 
represents product and ES is the intermediate enzyme- 
substrate complex. 

In the present model, the reaction is limited by diffu- 
sion only and the formation of the transition state com- 
plex ES is not taken into account. The conversion of 
substrate to product is also assumed to be instantaneous. 
Note that diffusion has the slowest time scale in this prob- 
lem. Hence, the above enzymatic reaction reduces to an 
irreversible one as 



E + S -> E + P. 



(2) 



The final equilibrium state corresponds to conversion of 
all substrates to products. A Monte Carlo (MC) al- 
gorithm has been developed to study diffusion limited 
enzymatic reaction as in Eq[2] in the presence of inert 
macromolecules. The algorithm is developed on the 2d 
square lattice of size L x L. For simplicity, the motion of 
the macromolecules is ignored and these act as immobile 
and inert obstacles. The inert obstacles do not interact 
with cither among themselves or with the substrate or 
product. The obstacles (O), enzyme (E), substrate (S) 
and product (P) are all represented as point particles 
in this model. It is also assumed that there exists only 
one immobile enzyme in the whole system. The enzyme 
is placed at the center of the lattice. After placing the 
enzyme, the obstacles and the substrates are distributed 
randomly over the lattice sites with their specified con- 
centrations Co and Cs respectively. A random number 
r is called from a uniform distribution of random num- 
bers between and 1 corresponding to each lattice site. 
If r < Cs, the site is occupied with a substrate and if 
Cs < r < a,f the site is occupied with an obstacle where 
cif = Cs + Co is the area fraction. The excluded volume 
condition is maintained, i.e., at any instant of time one 
lattice site cannot be occupied by more than one molecule 
of the same or different species. The substrate molecules 
diffuse through the space not occupied by the obstacles 
which will be referred as empty space later. As soon as a 
S reaches E, a product P is produced with unit probabil- 
ity. The diffusion of substrate or product in the system 
is modeled by simple random walk in presence of obsta- 
cles or disorder. At each MC time step, all the random 
walkers (all S and P) make an attempt to move to one of 
their nearest neighbors. The destination site, a site out 
of the four neighbors, of a random walker is chosen call- 
ing a random number r uniformly distributed between 
and 1. With respect to the present site, the destination 
site is going to be on the left if < r < 1/4, it is at the 



top if 1/4 < r < 1/2, it is on the right if 1/2 < r < 3/4, 
and it is at the bottom if 3/4 < r < 1. The destination 
site could be either empty or occupied by S, P, O or E. 
Depending on the status of the destination site, there are 
then four possibilities: (a) if the destination site is empty, 
the present S or P moves to the destination site, (b) if the 
destination site is occupied by a S or P, S or P remains 
on the same site, (c) if the destination site is occupied by 
an O, P or S also remains on the same site, and (d) if the 
destination site is occupied by the enzyme E, P remains 
on the same site but S is converted to P with unit prob- 
ability. If all the molecules of S and P are checked for an 
attempt of motion, time t (the MC time step) is increased 
to t+1. To ensure percolation of the substrate molecules, 
the maximum area fraction ay = Cs + Co is taken as 0.4, 
far below the percolation threshold. Note that, the per- 
colation threshold on the square lattice is w 0.59 [12j. 
Note that, the present non-equilibrium diffusion limited 
enzymatic reaction model is substantially different from 
that of lattice gas model incorporating equilibrium reac- 
tion rates proposed by Schnell and Turner 0] which leads 
to an unusual equilibrium constant equal to zero in the 
crowed environment [14j |. 

Cyclic boundary condition has been applied in the mo- 
tion of S and P. The simulation has been performed upto 
10 6 MC time steps on a 256 x 256 square lattice. The data 
are averaged over 100 ensembles. The time evolution of 
the system morphology for a/ = 0.1 with Cs = 0.01 
is shown in FigQ]at three different time. The black dots 
represent the substrates and the gray boxes represent the 
products. For clarity obstacles are not shown. It can 
be seen that the initial black dots are converted to gray 
boxes at the end. That means, the substrate molecules 
are diffusing, reacting with the enzyme, and are getting 
converted into products. In time, almost all the sub- 
strate molecules are converted to products and the prod- 
uct molecules also diffuse and spread all over the space 
uniformly. Lin and coworkers [15j simulated some ele- 
mentary kinetic reactions like A + B — ► with no ob- 
stacles under reflective boundary condition and observed 
Zeldovich crossover (segregation of A and B) Such 
segregation is not observed with periodic boundary con- 
dition in the present simulation. Effect of impenetrable 
boundary on diffusion limited reaction like A + A — > 
leads to different behavior depending on different bound- 
ary conditions 17 1 . 



RESULTS AND DISCUSSION 

Classical diffusion of a tracer particle in disordered sys- 
tems has already been studied extensively and the results 
are well understood [3] • Generally the diffusion is mod- 
eled by random walk and the disordered system is mod- 
eled by spanning percolation clusters [l2|. For studying 
diffusion, a quantity of interest is the root mean square 
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FIG. 1: The system morphology on a 256 x 256 square lattice is shown at three different times (a) t — 2 12 , (b) t — 2 18 and (c) 
t = 2 20 for substrate concentration C s = 0.01 and area fraction a/ = Cs + Co =0.1 (Co = 0.09). The black dots represent 
the substrates and the gray boxes represent the products. For clarity obstacles are not shown. The enzyme is represented by a 
cross at the center of the lattice. Products are formed due to the enzymatic reaction and in the long time limit almost all the 
substrates are converted into products. 



(rms) distance r(t) covered by the diffusing particle in 
time t. The rms distance r(t) in 2d is given by 



2/,- 



(3) 



where T> is the diffusivity of the system. The exponent 
k has a value 1/2 for diffusion on a regular lattice in the 
t — > oo limit. On the percolation cluster, diffusion is 
found to be anomalous and the value of k becomes less 
than l/2jlij|. The enzyme kinetic reaction inside a cell 
cytoplasm involves (i) diffusion of a large number of sub- 
strate molecules through the random structure of inert 
macromolecules, (ii) reaction with the enzyme to have 
products, and (Hi) finally diffusion of products from the 
enzyme through the same macromolecular crowding. The 
diffusion process involved here is then a collective mo- 
tion of a large number of particles in presence of disorder 
which is a complicated process than diffusion of a single 
tracer particle in a disordered medium. Self-diffusion is 
expected to play a nontrivial role along with the diffu- 
sion of S or P in presence of disorder in these systems. 
In order to check whether the enzyme kinetic reaction 
considered here is diffusion limited or not, one needs to 
analyze the the diffusive behavior of either the substrates 
or the products. To calculate the average diffusion length 
of the product particles, the coordinates {xi(t) 7 yi(t)} of 
each product i is recorded with time t. Time is measured 
starting from the birth of a product. The rms distance 
r(t) traveled in time t is then calculated as 
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FIG. 2: (a) Plot of diffusion length r(t) against time t in 
double logarithmic scale for different are fractions aj keeping 
substrate concentration at Cs = 0.01. (6) and (c) Plot of the 
local exponent k t versus time t for Cs = 0.1 and Cs ~ 0.01 
respectively. The same symbol set of (a) has been used in (6) 
and (c) for different area fractions a/ = Cs + Co- 



where (xo,yo) is the coordinate of the enzyme at the 
center of the lattice, N p (t) is number of products of age 
t. The data is then sample averaged over 100 ensembles. 

In FigJ^a), r(t) is plotted against time t in double log- 
arithmic scale for different area fractions df = Cs + Co 



keeping the substrate concentration constant at Cs — 
0.01. It can be seen that the magnitude of the diffusion 
length decreases with increasing a/. The collective mo- 
tion of the particles is then affected more and more by 
the presence of increasing inert macromolecules in the 



4 



system. However, in order to check the diffusive behav- 
ior of the particles one needs to estimate the exponent k 
defined in EqfSJ The local slope k t = d\og 2 r(t)/d\og 2 t 
of the curve log 2 r(t) versus log 2 1 is determined by em- 
ploying central difference method. In Fig[^6) and (c), k t 
is plotted against time t for two different substrate con- 
centrations Cs = 0.10 (6) and Cs = 0.01 (c) for the same 
set of area fractions a/ as in Figl^a). The value of k t 
saturates to 1/2 starting from a smaller value as t tends 
to a large value. Thus, a crossover from sub-diffusive to 
diffusive behavior has occurred for all area fractions in 
the long time limit. In the case of low substrate con- 
centration C s = 0.01 and high area fraction a/ = 0.4, kt 
shows certain anomalous behavior. Note that, at this pa- 
rameter regime the macromolecular concentration is 0.39 
which is just below 1 — p c ss 0.41 since the percolation 
threshold for a 2d square lattice is p c s» 0.59. The empty 
sites provides the connectivity for the substrate molecules 
all over the lattice. However, p c is defined on a infinitely 
large system. For a smaller system, even at the concen- 
tration of 0.39 the connectivity of empty sites may be lost 
for some of the ensembles considered. Consequently, the 
product may be trapped in a localized region around the 
enzyme and this may be the reason behind the anomalous 
behavior observed in this parameter regime. 

Since the parameter regime here is limited by diffu- 
sion, the enzyme kinetic reaction is then expected to be 
diffusion limited. Due to the enzyme kinetic reaction 
(given in Eq[5J) the substrates are converted to prod- 
ucts in time with unit probability on their encounter. 
In order to characterize the enzyme kinetic reaction, the 
number of products Np are counted as function of time 
t, the MC time step, for different substrate concentra- 
tions Cs and area fractions a/ = Cs + Co- In Fig|3] 
the product numbers Np is plotted against time t for 
different area fractions a/ with Cs = 0.01. Initially, Np 
increases linearly, then slows down and finally saturates 
in the long time limit. For low area fraction, it can be 
seen that the reaction is almost complete i.e.; most of the 
substrates given initially, Ns(0) — Cs X 1? ss 655, are 
converted to products exponentially as in classical equi- 
librium Michaelis-Menten kinetics though in the present 
model a non-equilibrium kinetics is considered. However, 
note that there is a considerable decrease in the product 
yield with increase in area fraction and their profiles are 
found not to follow an exponential increase. It has al- 
ready been predicted by numerical simulations that clas- 
sical Michaelis-Menten kinetics may not apply to enzy- 
matic reactions in crowded media [191]. In a Id model of 
reaction diffusion with disorder, Doussal and Monthus 
[20I ] also found large time decay in the species density 
via real space renormalization group calculations. The 
macromolecular crowding then could have a considerable 
and nontrivial effect on the enzymatic reaction rate. 

Initial rate of enzymatic reactions determines most of 
the molecular process and thus is an important quantity 
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FIG. 3: Plot of number of products Np versus time t for 
different area fractions 0/ = Cs + Co keeping substrate con- 
centration constant at Cs = 0.01. 



to estimate. Since non-equilibrium enzymatic reaction 
is considered here, the reaction rate R is defined as the 
ratio of the number of products Np to time t for 10% 
conversion of the substrates. R is then sample averaged. 
A similar analysis has also been performed for Np ver- 
sus t plots corresponding to Cs = 0.1 for different area 
fractions aj. In FigHJa), the normalized reaction rate 
R n = R/Cs is plotted against obstacle concentration Co 
for two different substrate concentrations Cs = 0.01 (cir- 
cles) and Cs = 0.1 (squares). Note that, area fraction 
aj = Cs + Co is not a good parameter to study the re- 
action rate since a/ will remain finite for finite Cs even 
at Co = 0. In the inset, R n is also plotted against Co 
in semi logarithmic scale. There are few things to notice. 
First, the reaction rate is decreasing with the increase in 
obstacle concentration Co in a nonlinear fashion. Sec- 
ond, the reaction rates are different for a particular Co 
even after normalizing by the substrate concentration Cs- 
Third, there is a monotonic decrease of ln(i?„) for small 
Co and deviates from linear decrease for large Co- The 
decrease in reaction rate with increasing crowding con- 
centration is expected and also observed in experiments 
[loL Eil ]. However, the dependence of the rate on the 
crowding concentration is different form the prediction 
made by Minton || in the transition state as well as dif- 
fusion limited enzymatic reaction in which a hump in the 
ln(i?„) versus Co curve is expected for an intermediate 
Co- Fourth, the normalized reaction rate is going to zero 
as Co approaches 1 — p c « 0.41. Beyond Co — 0.41, the 
obstacles could block the spanning clusters of the empty 
sites. Consequently the enzymatic reaction will be local- 
ized and the reaction rate is expected to go to zero. 

The above observations can qualitatively be under- 
stood in terms of diffusion and percolation phenomena. 
As Co increases, diffusivity is expected to decrease be- 
cause of the crowding due to obstacles. The influence of 
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FIG. 4: (a) Plot of normalized reaction rate R n — R/Cs 
against Co for two different Cs values 0.01 (circles) and 0.1 
(squares). ln(i? n ) is plotted against Co for the same Cs values 
in the inset. The same symbol set for different Cs values is 
used, (b) Plot of diffusivity D against Co for Cs = 0.01 and 
Cs = 0.1. The same symbol set of (a) is used. 



macromolecular crowding on the diffusion of solutes has 
been investigated in recent experiments utilizing different 
crowding agents and a reduced solute diffusion coefficient 
was observed with increa sing size and concentration of 
crowding macromolecules [2l| . An estimate of diffusivity 

V = (dr 2 (t) I dt) I 4 (as given in EqfS]) has been made uti- 
lizing the data of diffusion length r(t) for different sets 
of substrate (Cs) and obstacle (Co) concentrations. In 
FigHfb), V is plotted against Co for Cs = 0.01 (circles) 
and Cs = 0.10 (squares). Like reaction rate, diffusivity 

V is also studied as a function of obstacle concentration 
Co instead of a/. It can be seen that V also decreases 
with Co in a nonlinear fashion. First of all, it is inter- 
esting to note that the whole dependence of R n on Co 
is is in accordance with the behavior of V with Co ■ The 
enzymatic reaction rate in this parameter regime is there- 
fore mostly governed by diffusion and can be considered 
a purely diffusion limited enzymatic reaction. It is im- 
portant now to consider the low Co values, especially 
the case of Co = 0. For low Co values, T> is slightly 
less for Cs = 0.1 than that of Cs = 0.01 for the same 
Co • This slight decrease in V is due to diffusion through 
the self crowding at higher Cs- On the other hand, the 



reaction rate at zero obstacle concentration is expected 
to be proportional to Cs and V and R can be obtained 
as R ps Cs x T>. It can be seen that the normalized 
reaction rate R n obtained here is very close to the corre- 
sponding values of V at Co = for both the Cs values. 
At Co = 0, the self diffusion of the substrate molecules 
eventually determines the reaction rates and might be 
responsible for a slight decrease in R n for Cs = 0.1 with 
respect to Cs = 0.01 as seen in FigH^a). The effect of 
Cs in absence of obstacles has been verified numerically 
for several higher values of Cs and a considerable effect 
of self-crowding has been observed on the reaction rate 
as well as on diffusivity. Note that, R n values are slightly 
greater than V for almost all values of Co as it can be 
seen by comparing Figf4^a) and (b). This might have 
happened firstly due to the fact that the initial yield oc- 
cur only from the locally available substrate molecules. 
The diffusion length of these substrate molecules are very 
less in comparison to the expected diffusion length. Sec- 
ondly, one should note that the initial reaction rate for 
a given Co has to be calculated keeping the substrate 
concentration Cs fixed. However, in the present model 
the substrate concentration is decreasing with time as the 
substrates are being converted into products. The effect 
will be predominant for low Cs and small system size. 
Consequently the rate determination will be erroneous 
in the t — ► limit due to low yield. Hence, extreme care 
has to be taken in determining the initial reaction rate. 
The enzymatic reaction considered here is the completely 
diffusion limited and the results obtained are explainable 
in terms of diffusion in disordered systems. It is therefore 
intriguing to note that such a simple model of enzymatic 
reaction based on diffusion and percolation phenomena 
only, is able to explain qualitatively the experimental ob- 
servations [lO L ^2 1 as well as results obtained in compli- 
cated models HjiiSBI*!- Hence, diffusion is observed to 
be playing the crucial role in determining the enzymatic 
reaction rates. 

It should be emphasized here that enzymatic reactions 
occur in 3-dimensional space in living systems whereas 
the simulation is performed in 2-dimensions here. The 
simulation results obtained here agree qualitatively with 
the experimental observations and it is expected that the 
features of the model will be retained in higher dimen- 
sions also. The main difficulty in 3d simulation is in par- 
allel updating of a large number of substrate and product 
molecules during time evolution through a large number 
of MC time steps. Time required for the full conversion 
of substrate to product increases exponentially with the 
number of molecules (N$ — Cs X L d ) which strongly de- 
pends on the dimensionality of space for a fixed substrate 
concentration. However, for quantitative comparison of 
the results obtained in simulation with that of experi- 
ments, the model must be extended to three dimensions. 

The biochemical events in the densely crowded mito- 
chondrial matrix, the site for TCA cycle and fatty acid 
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oxidation pathway are largely governed by large macro- 
molecules of various sizes. It is thus important to in- 
vestigate the influence of crowding as exerted by macro- 
molecules of different sizes. A decrease in reaction rate 
has been observed in experiments for varying obstacle 
sizes keeping the obstacle concentration constant 22]. It 
seems that the complex interaction between obstacles 
and the substrate is size dependent and might be gov- 
erning the enzymatic reaction rate. It is expected that 
the diffusion of substrates across large macromolecules 
might be slow due to the complex interactions with the 
obstacles. In the present model of enzymatic reaction, 
this complex interaction between obstacle and substrate 
may be incorporated by introducing a residence time for 
the substrate molecules at each encounter with the obsta- 
cle. A slowing down in the reaction rate with increasing 
residence time has been observed in the simulation in ac- 



the experimental observations. 



cordance with the experimental results [22j . The details 
will be reported elsewhere. 



SUMMARY 

The effect of macromolecular crowding on the enzy- 
matic reaction rates has been modeled by a MC algo- 
rithm based on diffusion and percolation phenomena. 
The substrates, products, obstacles and enzyme all are 
represented by point particles. A single immobile en- 
zyme is considered and placed at the center of the lat- 
tice. The obstacles and the substrates are distributed 
randomly with their specific concentrations following a 
uniform distribution of random numbers between and 
1. The obstacles remain immobile throughout the simu- 
lation. It is found that the reaction is solely diffusion lim- 
ited under these conditions. The diffusion of substrates 
and products are modeled by a collective random walk. 
The products form gradually and subsequently almost all 
the substrates are converted into products after a long 
time. The initial reaction rate has been estimated for 
different substrate and obstacle concentrations. The nor- 
malized reaction rate has a nonlinear dependence on the 
obstacle concentration and found slightly dependent on 
the substrate concentration. The dependence of reaction 
rate on the substrate as well as obstacle concentrations is 
then qualitatively understood with the help of diffusion 
and percolation theory. The results qualitatively explains 
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